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A Wheeler- De wit t quantum constraint operator for four-dimensional, non-perturbative 
Lorentzian vacuum quantum gravity is defined in the continuum. The regulated 
Wheeler-DeWitt constraint operator is densely defined, does not require any renor- 
malization and the final operator is anomaly-free and at least symmmetric. 
The technique introduced here can also be used to produce a couple of completely 
well-defined regulated operators including but not exhausting a) the Euclidean 
Wheeler-DeWitt operator, b) the generator of the Wick rotation transform that 
maps solutions to the Euclidean Hamiltonian constraint to solutions to the Lorentzian 
Hamiltonian constraint, c) length operators, d) Hamiltonian operators of the mat- 
ter sector and e) the generators of the asymptotic Poincare group including the 
q-( quantum ADM energy. 

Attempts at defining an operator which corresponds to the Hamiltonian constraint 
> of four-dimensional Lorentzian vacuum canonical gravity |l| have first been made within 
the framework of the ADM or metric variables (see, for instance, 0) This formulation of 
Sh ; the theory seemed hopelessly difficult because of the complicated algebraic nature of the 
■ Hamiltonian (or Wheeler-DeWitt) constraint. It was therefore thought to be mandatory 
to first cast the Hamiltonian constraint into polynomial form by finding better suited 
canonical variables. That this is indeed possible was demonstrated by Ashtekar 
There are two, a priori, problems with these Ashtekar connection variables for Lorentzian 
gravity : 1) they are complex valued and are therefore subject to algebraically highly 
complicated reality conditions, difficult to impose on the quantum level, which make sure 
that we are still dealing with real general relativity and 2) the Hamiltonian constraint 
is polynomial only after rescaling it by a non-polynomial function, namely a power of 
the square root of the determinant of the three-dimensional metric, that is, the original 
Wheeler-DeWitt constraint has actually been altered. 

A solution to problem 1) has been suggested in |4| (see also ||) : namely, one can 
define real Ashtekar variables |J which simplify the rescaled Hamiltonian constraint of 
Euclidean gravity and then construct a Wick rotation transform from the Euclidean to 
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the Lorentzian regime where the complex Ashtekar variables are needed to simplify the 
Lorentzian constraint. The advantage of these real- valued variables is that they allow for 
the construction of a mathematically rigorous kinematical framework by means of which 
constraint operators can be regulated in a non-ambiguous fashion. In particular, this 
framework has been successfully employed to arrive at the complete solution of the Gauss 
and diffeomorphism constraints |7| . The main problem left is then to rigorously and con- 
sistently define the Hamiltonian constraint. Given the developments in connection with 
the above mentioned Wick transform, in order to define then the Lorentzian Wheeler- 
DeWitt constraint operator it would be sufficient to define the Euclidean operator and 
the generator of the transform. 

However, concerning problem 2), apart from the non-appealing fact that the Wheeler- 
DeWitt constraint was altered, what is worse is that the new Hamiltonian constraint 
carries density weight two. On general grounds, such an operator would need to be renor- 
malized (thus introducing a length scale) which breaks diffeomorphism invariance. A 
solution to this problem was first suggested in || : the idea is to take the square root of 
the rescaled Euclidean Hamiltonian constraint. While this seems to be a required techni- 
cal step to do in order to preserve diffeomorphism invariance there remain problems with 
it that have to do with taking the square root of an infinite number of non-commuting, 
non-positive, not self-adjoint operators. 

In this letter we introduce a novel technique which shows that 

• The requirement that the Hamiltonian constraint should be polynomial can be dis- 
carded. One can quantize the original Wheeler-DeWitt constraint in its non-polynomial, 
unrescaled form in a satisfactory way. The resulting operator is surprisingly simple and 
the problem of computing its kernel is conceivable. 

• Since the Wheeler-DeWitt constraint carries density weight one, problem 2) disappears, 
we never have to take ill-defined square roots. 

• We never introduce complex variables and we never need to deal with Euclidean gravity. 
A Wick rotation transform is completely unnecessary : problem 1) also disappears. Quite 
surprisingly, however, the unrescaled Euclidean Hamiltonian constraint operator and the 
generator of the Wick transform appear very naturally in our regularization procedure so 
that we cannot avoid to construct them simultanously as a side result. This could turn 
out to be important for the task of actually solving the Wheeler-DeWitt constraint. 

• It should be stressed at this point that all versions of Hamiltonian constraint operators 
for canonical quantum gravity in the continuum that have been constructed so far in 
the literature, whether in the connection representation || or in the loop representation 
|8], |TUJ, suffer from the fact that they provide quantizations only of the Euclidean theory. 
Their validity relies heavily on the assumption that one can actually exponentiate the 
generator of the Wick transform as outlined in 0] . The mathematical difficulties associ- 
ated with this step are highly non-trivial and it is even possible that the Wick transform 
simply does not exist. It is therefore the more important that we are able to work directly 
with the Lorentzian constraint for which the correct reality conditions are implemented 
right from the beginning. 

• The operator we obtain is not to be confused with the operator defined in [|TTJ (which 
is also formulated in terms of real variables) because a) we work in the continuum rather 
than on a lattice and b) our operator is completely well-defined on the whole Hilbert 
space while the one in All] suffers from singularities on a large subspace of it. 
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Let us fix the notation. Let the triad on the spacelike hypersurface £ be denoted by e l a , 
where a,b,c, ... are tensorial and i, j, k, ... are SU(2) indices. The relation with the intrinsic 
metric is given by q ab = e l a e{5ij. It follows that det(g) := &et((q ab )) = [det((e^))] 2 > 0. 
The densitized triad is then defined by Ef := det((e£))e" where e" is the inverse of e l a . We 
also need the field K l a = eK ab e b , e = sgn(det((e^))), where K ab is the extrinsic curvature 
of S. It turns out that the pair {K l a , E^) is a canonical one, that is, these variables obey 
canonical brackets {K l a (x) , E b (y)} = 5 l ^\x,y)8 b a 8 l j . 

Let the spin-connection (which annihilates the triad) be denoted by T l a . Then one can 
show that (A l a := r„ + K^E" 1 ) is a canonical pair on the phase space of Lorentzian 
gravity subject to the SU(2) Gauss constraint, the diffeomorphism constraint and the 
Wheeler-DeWitt constraint (neglecting a term proportional to the Gauss constraint) 
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H := Jdet(q)[K ab K ab - -{K a a f - R] = tr((F ab - 2R ab )[E a , E }) (0.1) 

where F ab and R ab respectively are the curvatures of the SU (2) connection A\ and the 
triad e l a respectively. 

What has been gained by reformulating canonical gravity as a dynamical theory of SU(2) 
connections is the following : if, as we do in the sequel, one makes the assumption that 
there exists a phase for quantum gravity in which the excitations of the gravitational 
fields can be probed by loops rather than, say, test functions of rapid decrease, then 
one has access to a powerful calculus on the space of (generalized) connections mod- 
ulo gauge transformations A/Q and, in particular, there is a natural choice of a diffeo- 
morphism invariant, faithful measure /io thereon which equips us with a Hilbert space 
7i := L 2 (A/Q,dfi ), appropriate for a representation in which A is diagonal. Moreover, 
Gauss and diffeomorphism constraints can be solved (see M and references therein for an 
introduction into these concepts). 

The remaining step then is to give a rigorously defined quantum operator corresponding 
to the Wheeler-DeWitt constraint and to project the scalar product on its kernel. We do 
this in a series of three steps. 

Step A) We begin by giving meaning to an operator corresponding to the Euclidean 
Hamiltonian constraint 

H E := -^=ti(F ab [E\ E b }) . (0.2) 
'det(g) 



The method applied in [§, |9|] is to absorb the prefactor 1/ y'det(g) into the lapse function 
and to give meaning to the operator corresponding to the square root of the trace. There 
is also an approach fl(| that avoids taking the square root, however, then one discovers 



a singularity which needs to be renormalized and this breaks diffeomorphism invariance. 
It can be recovered upon removing the regulator but then, to the best of our knowledge, 
all manipulations only have a quite formal character. Finally, in [JTTJ the Wheeler-DeWitt 



constraint is multiplied by a power of y det(g) to render it polynomial but such a procedure 
only works on the lattice where the density weight does not matter. In a final continuum 



limit one will ultimately encounter singularities of even worse character than in |10|. 



By employing the method described below we can avoid the complications that arise in 



all of these approaches. The presentation will be brief, details will appear elsewhere |L2 



We wish to impose the (self-adjoint) constraint operator on diffeomorphism invariant 
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distributions ip on A/G, that is, we evaluate ip on the action of the constraint on so- 
called "cylindrical" gauge invariant functions / [0] and require the resulting number to 
vanish, that is ip(Hf) = where t/>(/) = f^jgdfio^f- in brief terms, gauge invariant 
cylindrical functions on the space of (generalized) SU(2) connections are just finite linear 
combinations of traces of the holonomy around closed loops in S. Each such function 
therefore may equally well be labeled by the closed graph 7 consisting of the union of all 
loops involved in that linear combination. Such a graph consists of a finite number of 
edges ex, ..,e n and vertices v%, ..,v m . So, a function cylindrical with respect to a graph 7 
typically looks like f(A) = f y (h ei (A), .., h en (A)), where h e (A) is the holonomy along e for 
the connection A and / 7 is a gauge invariant function on SU(2) n . 
We are now ready to explain the main idea. 
• The total volume of S is given by 

V := f d 3 x^|det(g)| . (0.3) 



Since we will be using only the variation of V it is understood that if E is not compact 
then we first take a one parameter family of bounded subsets £# C X = where 
Tir C for R < R! to compute the variation of J- Sr \J\ det(g) | and then take the limit 
R — > 00. 

The first fact that we are going to use is that there is a well-defined, self-adjoint operator 



V on 7i corresponding to V [13] whose action on cylindrical functions is perfectly finite : 



(following the second reference in 




o 01 Yl e[ei,ej,eK)e ijk X\X 3 jX k K \ / 7 (#i,..,# n ) (0.4) 

° ' e / ne J ne Jf =D / 

where e(e/, ej, ex) = sgn(det(e/(0), ej(0), ex(0))) and ^(7) is the set of vertices of 7. We 
have abbreviated gi = h ei (A) and Xj = X(gj) is the right invariant vector field on SU(2) 
(we have chosen orientations such that all edges are outgoing at v). This demonstrates 
that V is a finite and well-defined operator on cylindrical functions. 
• Secondly, we exploit the elementary chain of identities (e abc has density weight one) 

[E a , E % _ e abc ee i^ = 2e abc JV = 2e abc^ yy _ ^ 



det(g) '" " SEfa) 

Then the Euclidean Hamiltonian constraint functional can be written (N is the lapse 
function) 

H E [N] = 2 f d 3 xN(x)e abc tr(F ab {A c , V}) . (0.6) 

We now triangulate E into elementary tetrahedra A and for each A we pick one of its 
vertices and call it v(A). Let ej(A), i — 1,2,3 be the three edges of A meeting at v(A). 
Let otij(A) := Ci(A) o a^- (A) o ej(A) -1 be the loop based at v (A) where is the obvious 
other edge of A connecting those endpoints of e», e^- which are distinct from f (A). Then 
it is easy to see that 

H£[N] := -^^^(^(^^^{/i- 1 ^, V}) (0.7) 
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tends to 2 J A Ntr(F A {A, V}) as we shrink A to the point v(A), N v := N(v(A)). More- 
over, H^[N] is gauge- invariant. 
Let the triangulation be denoted by T. Then 

H*[N] = H%[N] (0.8) 



is an expression which has the correct limit (|0.6| ) as all A shrink to their basepoints (of 
course the number of tetrahedra filling S grows under this process). 
The reason for doing this is clear : If we now simply replace V by V and the Poisson 
bracket by 1/ih times the commutator then 

9/V 

H*[N] := £ Hi[N], Hi[N] := -^^{^^[h^, V}) (0.9) 

is a regularized operator with the correct classical limit and whose action on cylindrical 
functions is indeed finite ! Namely, as may be suspected from the expression ( p.4j ), we 
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H§[N]f= £ Hi[N}f„ (0.10) 

AeT; Ari7^0 

that is, a tetrahedron contributes to the action on / only if it intersects the graph. More- 
over, as one can show in more detail, it only contributes if it intersects the graph in one 
of its vertices. Therefore, if we choose our triangulations T such that for each point of S, 
the number of edges of tetrahedra of T intersecting it, is uniformly bounded (uniformly 
in T) from above by some integer then, no matter how fine the triangulations is, there 
are always only a finite number of terms involved in the sum ( |0.10| ). 



It is amazing that one got expression (p.9|) almost for free once one knows that the volume 
operator is well-defined on holonomies, no ill-defined products of distributions arise, we 
do not encounter any singularities, no renormalization of the operator is necessary. 
One can determine, for each graph, a triangulation which is diffeomorphism covariant |i2| , 



meaning that the prescription of how to attach the loops ay(A) moves with the graph 
under diffeomorphisms in an appropriate sense. Then, for reasons first observed in ||, 
when one evaluates a diffeomorphism invariant distribution ip on ( |0.10| ), the number one 



gets depends only on the diffeomorphism class of the loop assignment aij(A). Therefore, 
in this diffeomorphism invariant context, the loops a^A) can be chosen as "small" and 
the triangulation as "fine" as we wish, the value of ( p.lOj ) on ip remains invariant and in 
that sense the continuum limit has already been taken. 

However, the operator ( p.9| ) still carries a sign of the regularization procedure : it depends 
on the diffeomorphism class [T] of the triangulation assignment which labels the freedom 
that we have in our regularization scheme. It is therefore not an entirely trivial task 
to check whether our operator H^[N] is anomaly-free, meaning that [H^[M], H^[N]]f 
vanishes for any cylindrical function / and lapse functions M, N when evaluated on a 
diffeomorphism invariant state. This severely constrains the freedom that we have in our 
choice of [T] . Specifically, a solution to the anomaly freedom problem is obtained if all 
the loops otij(A) are chosen to be kinks with vertex at v(A) ! That is, the arc Oij(A) joins 
the endpoints of e*(A), e 3 -(A) in an at least C l fashion. An arbitrary attachment of ciy(A) 
is insufficient to guarantee anomaly-freeness. In order to prove that the commutator van- 
ishes we need to use diffeomorphism invariance as follows : notice that the action of the 
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Euclidean Hamiltonian constraint is actually defined only up to a diffeomorphism. We 
need to make sure that for each choice of loop assignment within the same diffeomorphism 
class for either of the operators H E (M), H E (N) and their products the commutator van- 
ishes when evaluated on a diffeomorphism invariant state ip. When one performs the 
calculation it turns out that for each such choices one gets a sum of expressions of the 
form k(M, N)[U(<p)f — U(<fi')f] where k(M,N) is a certain function depending on the 
lapses M, N only, /' is a function cylindrical with respect to a graph which is bigger than 
the one that / depended on, 0, <fi' are certain diffeomorphisms depending on our choices 
of loop assignment and finally [/(</>) / 7 = f<j,M is a unitary representation of the diffeomor- 
phism group Diff(E) on TC. It is now obvious that ip vanishes on each of these expressions 
separately. 

Finally, upon taking a symmetrical ordering of ( P-9| ) we manage to arrive at a symmet- 
ric operator|]. Note that we have no problems in ordering V to the left or to the right 
of the holonomies involved since V has a finite action on holonomies of A as is clear 
from (|0.4| ) . The constraint algebra remains non-anomalous even after symmetric ordering 
which seems to be in conflict with general arguments raised in flTSl for finite-dimensional 



models and which show that in symmetric ordering the constraint algebra does never close 
with the generator of the diffeomorphism group appearing to the right of the structure 
functions. The resolution of the apparent contradiction is related to the fact that the 
infinitesimal generator of U((fit), where (fit is a one parameter subgroup of Diff(S), cannot 
be defined since U((f>t) does not act strongly continuously on H [[F]. Therefore the question 
of whether the diffeomorphism constraint operator appears to the right in the expression 
of the commutator of two Hamiltonian constraints cannot even be asked. 

Step B) Recall that the integrated (densitized) trace of the extrinsic curvature is up to 
a constant factor just the time derivative of the total volume with respect to the integrated 
Hamiltonian constraint (which is a signature invariant statement) : 

K := J d 3 x^/det{q)K ab q ab = -J d 3 xK l a E? = -{V, H E [N = 1]} (0.11) 

which also can be readily verified. Then we are naturally led to define 

k T :=^\V,H*[l]]. (0.12) 

If we order H E appropriately then both operators on the right hand side are symmetric 
and finite and we have produced a finite and symmetric expression for Kt- 

Step C) The aim is to write down an operator version of the Wheeler-DeWitt constraint 
(p.l|) which is perfectly well-defined and finite. We have two strategies at our disposal. 

Strategy 1) 
Upon realizing the following identity 

SK 

Kl = — = {Al,K}, (0.13) 



1 actua lly a little more care is needed here Jl^] : we need to make sure that the family of operators 
given by ( |3.10 ) is consistent |Q, that is, they all are projections to cylindrical subspaces of Ji of one and 



the same operator on Ti 
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the developments around step B) motivate to get rid of the complicated curvature term 
R a b involved in (p.l|) in favour of K. We have 

H + H E = -JL=tr([K a ,K b ][E a ,E b ]) = —=J=tr([{A a , K}, {A b , K}} [E a , E b ]) 



det(g) v /det(g) 
= 4e abc ti([{A a , K}, {A b , K}]{A C , V}) = 8e abc tr({A a , K}{A b , K}{A C , V}) . (0.14) 

The last identity suggests to define the regularized Wheeler-DeWitt operator in complete 
analogy with the Euclidean Hamiltonian constraint operator in the following, manifestly 
gauge invariant way 



H T [N] :-- 



(0-15) 

and we see that the Wheeler-DeWitt constraint operator can be built alone from the 
volume operator V and operators corresponding to holonomies along the edges of the 
tetrahedra of a triangulation. 

As it stands it is not self-adjoint yet but it is clear that a symmetrical ordering can be 
performed (without picking up singularities) to render it symmetric and we expect it to 
possess self-adjoint extensions. 

Finally, for the same reason that H E [N] is anomaly- free, Ht is anomaly free as well. 
Strategy 2) 

This strategy simplifies the problem of finding solutions to the constraints as ( p.9| ) is less 
complicated than ( p,15| ). It is a luxury at our disposal which we may use or not, however, 
it is not a necessary step. 

The generator of the Wick rotation transform can be defined now as Ct '■= (tt/2)Kt and 
as argued in f|] we may just define 

H T ■= Wt 1 HtW T) where W T := exp(-l/hC T ) . (0.16) 

We would proceed by first finding solutions to H E ip E = and then just analytically 
continue them to find solutions to the Lorentzian Wheeler-DeWitt constraint in a 
holomorphic representation. Such solutions are mapped unitarily [|J to solutions ip := 
Ut tp<c = (Wt^^Pe to H T ip = in the real representation. Expression (P.16|) has the 



disadvantage that whenever is symmetric, Ht is not. The motivation for having Ht 
symmetric is because we wish to find its kernel in the form of generalized eigenvectors . 
Expressions ( p,15|) , ( p,16|) are anomaly-free, densely defined operators corresponding to 



the original Wheeler-DeWitt operator in the continuum for non-perturbative four-dimensional 
Lorentzian canonical quantum vacuum gravity. 

Some final comments are in order : 
• Not even our Euclidean operator (|0.9|) and the ones proposed in |I| and || || have 
anything to do with each other, they are entirely different. The only thing they share 
is that their classical limits or the square thereof are proportional to each other. It is 
therefore to be expected that the Euclidean solutions that have been found already in 



the literature for the operators defined in [10, Bl are far from being annihilated by our 



operator. What speaks for our operator is that a) none of the operators in f|, 0] can 
be used to define K along the lines proposed here because it was crucial that the classical 
identity K = —{V,H E [1]} holds and b) our procedure leads to the quantization of the 
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original Wheeler-DeWitt constraint, rather than a modified version thereof. 

• There is a lot of freedom involved in the regularization step reflecting the fact that the 
quantum theory of a given classical field theory is not unique. An important, unresolved 
question is how to select the correct (or physically relevant) regularization of H . 

• As the Hamiltonian constraint operator on a given graph 7 reduces to a finite number 
of mutually commuting (in the diffeomorphism invariant context), symmetric constraint 
operators (which we expect to have self-adjoint extensions), one for each vertex of 7, 
we can exponentiate it and it seems feasable to determine the space of diffeomorphism 
invariant solutions to the Wheeler-DeWitt constraint as well as a physical inner product 
thereon by the group averaging method |16|, [7]. To illustrate this, let the function / 7 be 
cylindrical with respect to a graph 7 and denote by ^(7) its set of vertices. Then the 
constraint equation has the structure HT[N]f 1 = Y,vev(-y) N(v)HT, v f~/ = and is formally 
solved by (fi v is the Haar measure on the Abelian group generated by Ht,v) 

[/] 7 := J] / 'dfi v (a v )e ia ^ / 7 (0.17) 
vev(-y) 



Next, we take the group average |7j over the graphs on which the decomposition of ( p.17 ) 



into cylindrical functions depends to obtain a solution [f][ 7 ] to both, the diffeomorphism 
and the Hamiltonian constraint and the physical inner product would be 

< [f][7h lf%'] > P hys- = < [f}[ 7 ] J!y> > (0.18) 

where the second inner product is the one on 7i. The construction of interesting observ- 
ables would parallel related procedures displayed in J?J . 



• The final expression of the Wheeler-DeWitt constraint (|0.15|) is surprisingly simple : on 
each cylindrical function it is a low order polynomial in the volume operator and holon- 
omy operators and therefore one can hope to find exact solutions. While no solution could 
be found until now in closed form (except for cylindrical functions on two-valent graphs 
which, however, do not take the presence of the curvature term F a \, fully into account) 
the intuitive picture that arises concerning the action of the Wheeler-DeWitt constraint, 
is as follows : Recall [[17]] that a spin-network state is a gauge invariant function T^j(A) 



cylindrical with respect to a graph 7 where the dependence on its edges e 1 , ..,e n is through 
the matrix elements of irreducible representations of SU(2) labeled by spins j = ji, .., j n , 
evaluated at the holonomy along the corresponding edges. On such spin-network states 
the Hamiltonian constraint acts by annihilating, creating and re-routing the quanta of 
angular momentum associated with the edges of the graph in units of ±h,±h/2,0. Re- 
markably, the spectrum of the Hamiltonian constraint operator at a given vertex is largely 
determined by the spectrum of the volume operator so that it becomes of utmost impor- 
tance to gain control over it [JT8] . 



• It is clear that the method proposed here opens access to other well-defined and finite 
operators which were so far out of reach in a representation in which the intrinsic metric 
is not diagonal because we are able to make sense out of an operator corresponding to 
q a b. Examples are operators corresponding to the length of a curve [[Bl, matter Hamilto- 



nians for canonical Yang-Mills theory or the ADM Hamiltonian pi| . It is extremely 
interesting to see whether the latter Hamiltonian is at least positive semi- definite. 
• We define the spin-network representation to be the abstract representation defined by 
< A\j,j >:= T 7 j(v4) where as usual < A'\A >= 6 IM3 (A', A). Then it can be shown that 
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the ADM Hamiltonian acts effectively by multiplying the states > by certain alge- 
braic factors depending on the spins j. In other words, the spin-network representation 
can be interpreted as the "non-linear Fock-" (or occupation number) representation for 
quantum gravity Roughly, the spins associated with the edges of a graph indicate the 
"number" of elementary string-like excitations of the gravitational field along the edges 
of the graph very much like the integers indicate the number of photons associated with 



a certain array of excited modes for QED (see |]T2| for a more detailed discussion). 
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